Critical Swelling of Particle-Encapsulating Vesicles 
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We consider a ubiquitous scenario where a fluctuating, semipermeable vesicle is embedded in 
solution while enclosing a fixed number of solute particles. The swelling with increasing number 
of particles or decreasing concentration of the outer solution exhibits a continuous phase transition 
from a fluctuating state to the maximum-volume configuration, whereupon appreciable pressure 
difference and surface tension build up. This criticality is unique to particle-encapsulating vesicles, 
whose volume and inner pressure both fluctuate. It implies a universal swelling behavior of such 
vesicles as they approach their limiting volume and osmotic lysis. 

PACS numbers: 87.16.dj, 64.60.Cn, 64.60.an, 64.60.F- 



Membrane vesicles are fluctuating closed surfaces of 
0.1-10 /im scale, made of a flexible bilayer of amphiphilic 
molecules in aqueous solution. Serving as a simple model 
of biological compartments (e.g., red blood cells), they 
have been one of the most extensively studied soft-matter 
systems [H, H, H[. Numerous works have addressed the 
elasticity and statistical mechanics of the membrane un- 
der various constraints, such as area and enclosed volume 
(e.g., (3,0,0), or area and pressure difference across the 
membrane [a, 0, 0] , yielding various shapes and shape 
transformations. Actual vesicles are always immersed in 
solution and thus contain both solvent (water) and so- 
lute. Such biomolecule-encapsulating vesicles are ubiqui- 
tous in cell functions such as signaling and transport into 
(cndocytosis) and out of (exocytosis) the cell [9(. They 
are also used in various applications as microreactors or 
delivery vehicles for cosmetics and drugs (liposomes) [lj| • 

The hydrophobic core of the bilayer membrane hinders 
permeation of both water and solute molecules. Over suf- 
ficiently short time, therefore, the vesicle volume is fixed. 
Yet, while the activation barrier for water exchan ge i s 
of the order of 20T (T being the thermal energy) [llj . 
the barriers faced by the solute particles are typically 
much higher due to their size and/or charge, resulting 
in membrane permeabilities which are orders of magni- 



tude lower [11|. Moreover, water exchange across the 
membrane can be tremendously enhanced (indeed, bio- 
logically controlled) via aquaporin channels, which lower 
the barrier to below 8T |l2j. At sufficiently long times, 
therefore, most vesicles are found in a wide semiperme- 
able regime, where water can be considered as exchanged 
between the interior and exterior, while the solute re- 
mains trapped inside. As a result, it has been assumed 
that the exterior solution concentration and number of 
encapsulated particles determine the vesicle volume in 
practice [H, 0, [1| — the mean volume adjusts through 
water permeation so as to annul the osmotic pressure 
difference across the membrane. This scenario has been 
experimentally verified [T3| and utilized to measure mem- 
brane permeabilities of various solutes [ll| and create 
osmotic motors 14 1 . Volume fluctuations around the os- 



moticall y de termined mean value have been considered as 
well [la, Ufa- However, at high swelling, as the vesicle ap- 
proaches its maximum volume, this volume-adjustment 
description must break down, and appreciable pressure 
difference and surface tension will begin to build up. Fur- 
ther swelling eventually leads to pore formation and os- 
motic lysis [3, [I?], Eli ■ The change in swelling behavior, 
in particular, whether it is a smooth crossover or a sharp 
transition, is the focus of the current Letter. 

We describe the vesicle as a closed surface composed 
of N molecules and having maximum volume K max ~ 
a 3 TV 3 / 2 , a being a molecular length comparable to the 
membrane thickness. It is assumed that at V^ lax the 
vesicle has a nonextensive number of configurations. The 
vesicle encloses Q solute particles, which do not directly 
interact with the surface other than being trapped inside 
it. The vesicle is immersed in a solution of fixed concen- 
tration and temperature, which exerts an outer osmotic 
pressure p Q on the membrane. Since solvent molecules are 
exchanged between the interior and exterior, the vesicle 
volume is not specified and, hence, neither are its inner 
particle concentration and pressure. Thus, the partition 
function involves integration over all possible volumes, 



PoV/T 



Z(T, Po , Q,N) = J dVZ v {T, V, N)Z S (T, V, Q)e 

(1) 

where Z v and Z s are the canonical partition functions of 
the vesicle and solute particles, respectively, for a given 
volume V. The thermal energy T is hereafter set to unity. 
For the solute we write 



-QKQ/V) 



(2) 



where / is the canonical free energy per solute molecule. 

A key issue for the highly swollen vesicles studied here 
is how Z v behaves as V approaches V^ lax . It is shown 
below that, quite generally, 



Z V (V < Knax) ~ {V m 



V) 



olN 



(3) 



where a is a coefficient of order unity. This result read- 
ily follows from the two requirements, that (i) the vesicle 



2 



free energy be extensive in TV for V < V max , and (ii) the 
probability density function of volumes vanish at V maK . 
In more detail, Z v is found by integrating over all surface 
configurations the factor e~ H ^6(V — V[JZ\), where V[R] 
is the volume of configuration TZ, and H[TZ] its energy 
(including, e.g., contributions from bending rigidity and 
surface interactions). For V ~ VI nax one can generally 
represent the configurations by the amplitudes {u n } of A 
normal modes. (For example, in the case of nearly spher- 
ical vesicles these are spherical harmonics [H, fia. l20l|.) 
One then expands V[R] ~ V max — ^ n C n \u n \ 2 . Assum- 
ing that H is nonsingular at Vmax, and using the inte- 
gral representation of the delta function, we get Z v ~ 
e -H(v^) j d[u n ]d P e^[ip(V max -V-j:C n \u n \ 2 )}. Intc- 
gration over {u n } gives a factor of p~ x l 2 per mode which, 
upon integration over p, yields Eq. ^ with a = 1/2 [2lj |. 

Substituting Eqs. ([2]) and ([3]) in Eq. {I} while speci- 
fying the solute free energy /, one can perform the in- 
tegration over V for given Q, p Q , and N. In Fig. Q] we 
present the resulting mean volume, (V) = — dhiZ / dp Q , 
as a function of Q for an ideal solution, f(Q/V) = 
ln(Q/V) — 1. As A is increased, (V) is seen to approach 
a discontinuous first derivative at Q c — p V max . 
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FIG. 1: (Color online) Order parameter as a function of 
control parameter for an ideal solution encapsulated in vesi- 
cles of various sizes. Solid curves show the mean-field re- 
sults [Eq. ||SJ|], while dotted curves are obtained from numer- 
ical integration of the partition function [Eq. ft}]. Datasets 
from top to bottom (bottom to top in the inset) correspond 
to p = 1, N = 30 (green); p = 1, JV = 10 3 (blue); 
p = 5, JV = 10 3 (indigo); Po = 1, JV = 10 5 (red); and 
the JV — *• oo limit [Eq. (JBJ, black]. Only for the smallest vesi- 
cle size (JV = 30) are the mean-field and numerical results 
distinguishable. Inset shows rescaled data according to Eq. 
©, where the uppermost curve (solid black) is the theoretical 
scaling function. Values of p are in units of T/a 3 ; in all data 



sets a = 1/2 and V m 



i N 3/2 



lytically for a general (nonideal) solution. The parti- 
tion function can be rewritten as Z ~ f dVe~ F , with 
F = -aAln(^ nax - V) + Qf(Q/V)+ Po V. Minimizing 
F with respect to V and applying a first-order saddle- 
point approximation (which is equivalent to a mean-field 
assumption), we obtain the following equation for (V): 

Q 2 f(Q/(V))/(V)- Po (V) =aN(V)/(V max -(V)). (4) 



Expansion in V m: 
M = 1 - (V)/V m 



(V) yields for our order parameter, 



M 



N- 



y/[s-t(q)]*+4sg(q) - a - t(q)\ /(2g(q)) (5) 
^(\t\-t)/(2g), (6) 



where q = Q/V max is the solute concentration at V = 
Fm ax , t(p) = p 2 f'(p)/p — 1 is the rescaled difference 
between the solute pressure at concentration p and the 
outer pressure [t(p = q) acting as the control parameter 
of the transition] , s = aN/(p V max ) ~ JV" 1 / 2 , and<7(<7) = 
q 2 [2f(q)+qf"(q)}/p o >0 Q. 

Equations ([5]) and ^ describe a phase transition at 
q = q c which solves the equation t(q c ) = 0, i.e., for 
which, if the volume were equal to Vmax, the inner pres- 
sure would just balance the outer one. The parame- 
ter t(q) is related to the actual control parameter, Q 
or q, via the solute equation of state. In the ideal- 
solution example, f(Q/V) — \n(Q/V) — 1, the criti- 
cal point is at q c — p Q (Q c = p Vm ax ), and we have 
t = Q /Q c — 1 and g = 1. The transition occurs in the re- 
gion \t\ < A = (Asg) 1 / 2 ~ a- 3 / 2 P o 1/2 JV" 1 / 4 , along which 
M crosses over from finite values to very small ones, 



f-*/.9 - A^ ^! 1 <«-A 
M = I A/(2g) ~ N- 1 / 4 ^ |t| « A (7) 
[A 2 /(4gt) ~ JV- 1 / 2 *- 1 OA. 

In the thermodynamic limit, A — > oo, this crossover 
turns into a sharp corner [Eq. i.e., a discontinuity in 
dM/dt. From Eq. ([5]) we find that M follows a scaling 
law within the transition region, 



M/A = g^Mtf/A), M(x) = ( \J x 2 + 1 - x)/2, (8) 

which is verified in Fig. [T] (inset). In addition, we calcu- 
late from Eq. © the compressibility, \ — dM/dp a , 



We now proceed to investigate this criticality ana- 



X=(gp o r 1 x(t/A), x(x) = (1 - x/Vx* + l)/2. (9) 

Performing the next-order saddle-point calculation (i.e., 
including fluctuations beyond mean field) yields negligi- 
ble corrections to M, of order JV~ 3 / 2 , A" 5 / 4 , and A" 1 , 
respectively, in the three domains of Eq. ([7]). Thus, the 
mean-field description is accurate, as is also confirmed by 
numerical integration of the partition function (Fig. [I}. 
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Equation {H (upon division by (V)) is just the Laplace 
law, balancing the pressure difference across the mem- 
brane (left-hand side) with a surface term (right-hand 
side). We therefore identify the pressure difference and 
surface tension as 



Ap = (aN/V D 

RAp ~ aT 2 M 



-x 



a 



max)M- 

-2 H/T— l 



i- 3 N- 1/2 A-r\ 



(10) 



where R ~ aN 1 / 2 is the vesicle radius. Thus, Ap and 
a change from negligible values below the transition to 
appreciable ones above it. Specifically, a 3 Ap is of order 
Ar -i/2 j AT- 1 / 4 , and 1, while a 2 a ~ 1, TV 1 / 4 , and N 1 / 2 , 
below, at, and above the critical point, respectively. 

Note that Eqs. J]}-©, which underly the entire anal- 
ysis, contain no microscopic information. The model, 
therefore, is purely thermodynamic, in the sense that 
any specific microscopic model for the vesicle and en- 
capsulated solution (so long as the vesicle has a state 
of maximum volume and negligible entropy) should lead 
to the same results. (For example, inclusion of bending 
rigidity will merely change the prefactor in Eq. (J3J) . ) The 
invariance to the choice of model implies also that the 
continuous transition does not necessarily involve a di- 
vergent correlation length [Hj]. We have checked these 
statements for the specific example of a nearly spheri- 
cal envelope of N nodes and fixed total area 47T-R 2 ,, en- 
closing an ideal solution. The vesicle shape is defined 
in this case by R(9,p), the distance of the membrane 
from the center as a function of solid angle, whose de- 
viation from i?o can be decomposed into spherical har- 
monics, R{6,Lp) = Rq[1 + Y,l=^J2m=-l U lmYl m (8,^)], 

where (Z max + l) 2 = N. Integration of the resulting par- 
tition function over the amplitudes ui m within a saddle- 
point approximation recovers Eqs. (O-©. The correla- 
tion function, (uz m u;(- m )} ~ (M/N)/[l(l+l) — 2], I > 1, 
exhibits a critical suppression of amplitude but no di- 
vergent correlation length. Expectedly, this fluctuation 
spectrum is identical to that of a spherical membrane 
with surface tension a ~ Af _1 , in accord with Eq. (flC)]) . 

Equations 10) and (HJ) characterize the sharpening of 
the transition with increasing system size. If we recast 
them in the conventional finite-size scaling form [24| . 
M ~ RT 13 ^' M{R l ' v *t) and X ~ R l/v ' x(R 1/v " t), we 
readily extract /3 = 1, 7 = 0, and v* = 2. The values of 
(3 and 7 are consistent with the linear increase of M below 
the transition [Eq. ((6]) and Fig. [Tj . Notwithstanding the 
absence of a divergent correlation length, one can use v* 
to define a length scale, £ ~ a\t\~ l/ , such that the system 
lies in the critical domain if R < £. The divergent £ does 
not relate to correlations but to the competition between 
surface degrees of freedom (~ N) and three-dimensional 
ones (~ TV 3 / 2 ). This competition determines the width 
of the transition, A ~ [N/ (poKnax)] 1 / 2 , making it shrink 
with increasing N. Repeating the analysis for a ring in 
two dimensions yields a similar mean-field transition with 



identical exponents. We are not aware of another tran- 
sition whose mean-field limit has the exponents found 
above [25|. 

The phase transition just characterized is a unique 
feature of particle-encapsulating vesicles. It is a conse- 
quence of the effective inner pressure being dependent 
on the volume (through f(Q/V) for fixed Q), while the 
latter fluctuates. This leads to pressure difference and 
surface tension which are nonanalytic in Q [Eq. (|10p ] 
and a consequent breaking of the equivalence between 
the fixed-pressure (or fixed-tension) scenario and that 
of fixed Q. Indeed, if the enclosed solution is replaced 
with a given inner pressure p\ (i.e., upon substituting in 
Eq. ([I]) Z s = e PiV / T ), it is straightforward to show that 
M = aN[(pi - p )Knax] _1 , in agreement with Eq. (fTU|) 
[26l | . Therefore, in the case of a given pressure difference 
(or tension) the vesicle swells gradually with pi (or a) 
without criticality. Furthermore, replacing the particle- 
number constraint with a chemical potential [53] is equiv- 
alent (via the solute equation of state) to specifying the 
inner pressure. Hence, there is no criticality in the grand- 
canonical case either , an d the two ensembles are mani- 
festly not equivalent [281 ] . 

Another noteworthy limit is that of a pure solvent out- 
side the vesicle, p — > 0, where the current analysis yields 
Q c — > and A — * 00, i.e., the transition disappears. In 
this case the swelling of the vesicle toward its maximum 
volume occurs already for much smaller particle numbers 
Q, scaling with the area N rather than the volume (28l |. 

In summary, we have found that membrane vesicles, 
under rather general conditions, behave critically as the 
number of solute particles inside them is increased or, 
equivalently, the outer osmotic pressure is decreased. It 
should be possible to experimentally observe this phase 
transition, e.g., by creating vesicles and subsequently di- 
luting the outer solution in a controlled manner, or by 
using isotonic solutions of molecules with different mem- 
brane permeabilities [l8l |. We mention three points re- 
lating to such experiments. First, they should cover 
such time scales that the vesicle could be considered 
permeable to water. This can be sensitively controlled 
if water (aquaporin) channels are incorporated in the 
membrane, yet common lipid vesicles are also found in 
this re gim e over readily accessible time scales (~ 10 



13, 14 



18( | . Next we examine the assumption of 
a sharply defined maximum volume. One definition of 
Knax would be the volume enclosed by an unstretch- 
able vesicle of a given area once out-of-plane fluctua- 
tions have completely vanished 29]. This assumption 
should be relaxed when inplane (stretching) fluctuations 
become comparable to transverse ones. Since, for a tense 
membrane, the mean-square fluctuations of both modes 
have the same (quadratic) dependence on wavenumber, 
this crossover will occur simply when the surface tension, 
a ~ (T/a 2 )M -1 [Eq. (fTU)) ]. becomes comparable to the 
membrane stretching modulus, which is typically of or- 
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der 10 2 dyne/cm [3(J. For a ~ 1 nm this happens when 
M < 1CT 2 , i.e., when the mean volume deviates from 
Knax by less than 1%. Thus, we expect the transition 
from appreciable to small values of M, along with the cor- 
responding scaling behavior, to be manifest well before 
stretching becomes important. Since lipid membranes 
can sustain inplane strains of only a few percent be- 
fore rupturing [17[ , the crossover to stretching-dominated 
dynamics will be shortly followed by vesicle lysis [lij ]. 
Third, because of the weak dependence of the transition 
width on JV, A ~ (a 3 p /T)- 1 / 2 N- 1 ^, the observed be- 
havior will not be very sharp. A typical micron-sized 
vesicle has about N ~ 10 8 molecules in its membrane, 
leading, for a 0.1M solution, to A ~ 1CP 1 only. The crit- 
icality could be verified, nonetheless, by checking data 
collapse according to the scaling law, Eq. (JSJ). 

Thus, our assumptions concerning permeability, max- 
imum volume, and number of molecules do not rule out 
an experiment aimed at the predicted critical swelling. 
(The suppression of small fluctuations near the transi- 
tion, however, may be hard to resolve.) More generally, 
this study highlights the qualitative difference between 
semipermeable, particle-encapsulating vesicles and those 
having fixed volume or pressure. Since most natural and 
industrial vesicles belong to this class, their different be- 
havior should be taken into account, particularly in cases 
of high swelling and osmotic lysis. 
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